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ABSTRACT 

The Suzaku AGN Spin Survey is designed to determine the supermassive black 
hole spin in six nearby active galactic nuclei (AGN) via deep Suzaku stares, 
thereby giving us our first glimpse of the local black hole spin distribution. Here, 
we present an analysis of the first target to be studied under the auspices of this 
Key Project, the Seyfert galaxy NGC 3783. Despite complexity in the spectrum 
arising from a multi-component warm absorber, we detect and study relativistic 
reflection from the inner accretion disk. Assuming that the X-ray reflection is 
from the surface of a flat disk around a Kerr black hole, and that no X-ray 
reflection occurs within the general relativistic radius of marginal stability, we 
determine a lower limit on the black hole spin of a > 0.88 (99% confidence). 
We examine the robustness of this result to the assumption of the analysis, and 
present a brief discussion of spin-related selection biases that might affect flux- 
limited samples of AGN. 
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Introduction 



Ever since the seminal work of iPenrosd (Il969l ) and iBlandford fc Znajekl (jl977( ) , it has 
been recognized that black hole spin may be an important source of energy in astrophysics. 
Of particular note is the role that black hole spin may play in relativistic jets such as those 
seen in radio-loud active galactic nuclei (AGN) — the magnetic extraction of the rotational 
energy of a rapidly spinning black hole is the leading contender for the fundamental energy 
source of such jets. Indeed, it has been suggested that the spin of the central supermassive 
black hole (SMBH) is a crucial parameter in determining whether an AGN can form powerful 



jets (i.e., whether the source is radio-quiet or radio-loud: IWilson fc Co 



the accretion rate/mode must clearly have a role to play (jSikora et al. 



bert 



20071 ). 



1995|), although 



However, the importance of black hole spin goes beyond its role as a possible power 
source. The spin distribution of the SMB H population (and its dependence on SMBH 



mass) encodes the black hole growth history (iModerski fc Sikoralll996l : IVolonteri et al.ll2005 



Berti fc Volonterill2008l ). In essence, if local SMBHs have obtained most of their mass dur- 
ing prolonged prograde accretion events in a quasar phase of activity, or in major merg- 
ers with similar mass SMBHs, we would expect a population of rapidly rotating SMBHs 



(a > 0.6) due 1 


;o the 


RezzoUa et al. 


2008) 



the black hole and M is its mass. On the other hand, if mergers with much smaller SMBHs 



(iHughes fc Blandfordll2003[) or randomly-oriented accretion events of small packets of mate- 



rial ( iKing fc Pringld 120071 ) have been the dominant growth mechanism, most of the SMBHs 



would be spinning at a much more modest rate. 

To date, the cleanest probe of strong gravitational physics around SMBHs, including the 
effects of black hole spin, comes from examining relativistically-broadened spectral features 
that are emitted from the s urface layers of the inner accretion disk in response to irradiation 
by the hard X-ray source (IReynolds fc NowakI l2003l : iMillerl 120071). These spectral features 



have been observed and well-characterize d in both AGN (iTanaka et a. 



19951 ) and stellar-mass black hole systems (IMiller et al.ll2002 



Reis et al. 



..19951 : iFabian et al. 
2008|). The strongest 



feature in this so-called "reflection spectrum" is the fluorescent Fe Ka line (rest frame energy 
of 6.4 keV); in contrast to lines from other elements, its relative abundance, high energy and 
fluorescent yield make Fe Ka visible above the typical power-law continuum seen commonly 
in BH systems. Extreme Doppler effects and gravitational redshifts combine to give this line 



(Fabian et al. 


1989; 


Laor 


1991) 



XMM-Newton and Suzaku, combined w ith the latest models of reflection from an ionized 



accretion disk (e.g.. lRoss fc Fabianll2005l ) and variable-spin relativistic smearing models (e.g.. 
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Brenneman fc Reynolda l2006l : iDauser et al.l |2010| ) , are giving us our first glimpses at the 
spins of SMBHs. However, due to the high S/N required to characterize the subtle effects 
of SMBH spin, interesting spin con straints have only been determ ined for a small handful 
of AG N at present (11100-6 -30-15, iBrennernanfc Revnold^ bood; Fairall 9, Jschmoll et al. 



2009 



Gallo et al. 2011; Mrk 335 & NGC 7469 Patrick et al. 2010; see Table 2 



SWIFT J2 127.4+5654j Miniutti et al.lbood: 1H0 707-495. IZogh bi et al lboid Mrk 79 



Under the auspices of the Suzaku Key Projects program, we have initiated a series 
of deep quasi-continuous observations of bright, nearby AGN with the purpose of charac- 
terizing relativistic disk features in the spectra and setting constraints on the SMBH spin 
{Suzaku AGN Spin Survey; PI C. Reynolds). In this Paper., we present results from the first 
object to be studied under this program, the Seyfert 1.5 galaxy NGC 3783 {z = 0.00973; 
Theureau et al.lll998l ). This object possesses a high-column density and multi-component 
warm absorber that has been well-studied by every spectroscopic X-ray observatory, includ- 
ing a 900 ks campaign by Cfeandra using the High- Energy Transmissio n Grating Spectrometer 
fHETGS: E<aspi et al.l[2002[ IXrongold et~aDl2003[ iNetzer et al]l2003h . We show that, despite 
the presence of this complex warm absorber, refiection signatures from the inner accretion 
disk can be identified and characterized with sufficient accuracy to constrain SMBH spin. 
We conclude that the SMBH is rapidly spinning with a > 0.93 (90% confidence). This result 
is shown to be robust to the exclusion of the complex, soft region of the X-ray spectrum as 
well as to uncertainties in the XIS/PIN cross- normalization. 

This paper is organized as follows. Section |2] discusses the Suzaku observation of 
NGC 3783 and the basic reduction of the data. Section 3 then presents our modeling of the 
0.7 — 45keV time-averaged spectrum of NGC 3783, including our newly-derived constraints 
on the SMBH spin. Section 4 summarizes our conclusions on NGC 3783 and addresses the 
role of spin-dependent selection biases in AGN samples. 



2. Observations and Data Reduction 

NGC 3783 was observed by Suzaku quasi-continuously for the period 10-15 July 2009, 
with the source placed in the Hard X-ray Detector (HXD) nominal aimpoint. After eliminat- 
ing Earth occultations. South Atlantic Anomaly (SAA) passages, and other high background 
periods, the observation contains 210 ks of "good" on-source exposure. The XIS data (XIS 0, 
XIS 1 and XIS 3; XIS 2 has been inoperable since November 2006) were reprocessed using the 
xiapi script in accordance with the S«.at. ABC Guid^ along with the latent version of the 



^http://heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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CALDB (as of 29 March, 2010). XIS spectra and light curves were then produced according 
to the procedure outhned in the ABC Guide. For the XIS spectra, we combined data from the 
front-illuminated (FI) detectors XIS 0+3 data using the addascaspec script in order to in- 
crease S/N. The XIS spectra, responses and backgrounds were then rebinned to 512 spectral 
channels from the original 4096 in order to speed up spectral model fitting without compro- 
mising the resolution of the detectors. Finally, the XIS spectra were grouped to a minimum of 
25 counts per bin in order to facilitate robust fitting. The merged, background-subtracted, 
time-averaged FI spectrum has a net count rate of 4.960±0.002 cts s~^ for a total of 1.04 x 10^ 
counts. The total number of 2 — lOkeV counts is 6.26 x 10^. The total XIS 1 count rate is 
3.043 ± 0.004 cts s^i for a total of 6.40 x 10^ or 3.14 x 10^ when restricting to 2 - lOkeV. 
For all of the fitting presented in this paper, we allow for a global flux cross-normalization 
error between the FI and XIS 1 spectra. The XIS 1/FI cross-normalization is allowed to be 
a free parameter, and is found to be approximately 1.03. 

The HXD/PIN instrument detected NGC 3783, though the GSO did not. Data from 
PIN were again reduced as per the Suzaku ABC Guide. For background subtraction, we used 
the "tuned" non X-ray background (NXB) event file for July 2009 from the Suzaku CALDB, 
along with the appropriate response file and flat field file for epoch 5 data. The NXB 
background contributed a count rate of 0.2165 ±0.0004 cts s~^ to the total X-ray background 
from 16 — 45 keV. We modeled the cosmic X-ray ba ckground (CX B) contribution as per 



the ABC Guide, simulating its spectrum in XSPEC (lArnaudI Il996l ). The simulated CXB 
spectrum contributed a count rate of 0.0237 ± 0.0002 cts s^^ to the total X-ray background 
from 16 — 45 keV. The NXB and CXB files were combined to form a single PIN background 
spectrum. In comparison, the PIN data had a count rate of 0.4561 ± 0.0027 cts s~^ over the 
same energy range, roughly twice that of the total background. 

Because the PIN data only contain 256 spectral channels (vs. the 4096 channels in the 
unbinned XIS data), rebinning to 25 counts per bin was not necessary in order to facilitate 
fitting. Rather, we rebinned the PIN spectrum to have a S/N of 5 in each energy bin, 
which limited our energy range to 16 — 45keV. After reduction, filtering and background 
subtraction, the PIN spectrum had a net 16 — 45keV count rate of 0.360 ±0.002 cts s~^. We 
also added 3% systematic errors to the PIN data to account for the uncertainty in the non- 
X-ray background data supplied by the Suzaku calibration team. For most of the spectral 
fitting presented in this paper, we assume a PIN/XIS-FI cross-normalization factor of 1.18 
as per the Suzaku memo 2008-0«E However, for the Bnal Bts used to consttain the black 
hole spin in §3.4[ we investigate the effect of allowing this cross-normalization to be a free 



^http://heasarc. gsfc.nasa.gov/docs/suzaku/analysis/watchout. html 
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NGC 3783 




Time (s) 



Fig. 1. — Co-added and background subtracted XIS light curves in the soft (0.3 — 1 keV) and hard 
(2 — 10 keV) bands, together with the background subtracted PIN (16 — 45 keV) light curve. These 
light curves are shown with 5000 s bins. 

parameter. 

XIS light curves (both hard and soft band) as well as PIN light curves are shown in 
Fig. [H In the soft band (0.3 — 1 keV), the source is observed to undergo variability by a factor 
of almost two. Most of the large amplitude variability occurs on timescales of 50 — 100 ks, 
although there are occasional sharp fiares/dips that occur much more rapidly. Also note- 
worthy is that the amplitude of variability decreases as one considers higher-energy bands, 
suggesting "pivoting" of the spectrum about some energy above the Suzaku/PIN band. The 
detailed nature of this spectral variability will be the subject of another publication (Reis et 
al. , in preparation). For the remainder of this paper, we examine the high S/N spectrum 
from the time-averaged dataset. We restrict our energy range to 0.7 — lOkeV in the XIS 
data, ignoring energies below 0.7keV and from 1.5 — 2.5 keV to avoid areas of significant 
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(a) (b) (c) 




Energy (keV) Energy (keV) Energy (keV) 



Fig. 2. — Left panel: Simple power-law fit to the 3.5 — 45 keV XIS-FI+PIN spectrum. Middle panel: 
Zoom-in on the 4-8 keV region of the simple power-law fit. Note the probable "Compton shoulder" 
on the immediate low-energy side of the strong 6.4 keV emission line. Right panel: Residuals 
remaining when the broad iron line component is removed from a simple phenomenological fit to 
the 3.5 — 45 keV data (see §3.ip . In all panels the black points correspond to XIS 0+3 data, red to 
XIS 1 data and blue to PIN data. The green line represents a data-to-model ratio of unity. 

deviation between the three detectors, i.e., regions of known calibration uncertainty. 



3. Analysis of the time-averaged Suzaku spectrum 

3.1. A First Look at the Hard-band Spectrum 

It is instructive to begin by examining the hard-band (3.5 — 45 keV) XIS+PIN spectrum. 
A simple power-law fit to this band reveals significant spectral complexity (Fig. |2^,b). A 
narrow Ka fluorescence line of cold iron (6.4 keV) dominates; structure redward of this line 
indicates a possible Compton shoulder as well as an extended tail reaching down to ~ 4keV. 
Both the narrow iron line and the broad red wing likely originate from X-ray reflection 
and, hence, the convex spectrum between 8 — 40 keV is readily interpreted as the associated 
Compton reflection hump. Structure above the 6.4 keV line indicates a strong absorption 
feature at ~ 6.6 keV and/or an emission line at ~ 7keV (likely corresponding to a blend of 
the K/3 line of cold iron and the Lya line of Fexxvi). 

Guided by these identifications, we construct an heuristic model of the hard spectrum 
consisting of a power-law continuum, a narrow Fexxvi emission line (modeled as a Gaus- 
sian line centered at 6.97 keV with a = 10 eV), reflection from distant, low- velocity, cold 
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matter (described by the pexm on modej^) , and a relativistically broadened cold iron Ka line 



(described by the laor model; lLaorlll99ll ). This model produces an excellent fit to the data 
ix^/^ = 573/544(1.05)) with the following parameters: photon index T = 1.68tooi, refiec- 
tion fraction TZ = 0.87^q^I, emission line equivalent widths VFpeXxvi = 28;!^5eV, Wbroad = 
263^1 eV, inner edge of line emitting disk rin = 3.0^q § r^, index of line emissivity across disk 
q = 3.3ltQ 09, and disk inclination i < 9°. If we replace the pexmon model with a pexrav and 
three separate Gaussian lines for narrow Fe Ka (6.4 keV, a = 0.015 keV), Fe K/3 (7.06 keV, 
a = 0.015 keV) and Compton shoulder (6.25 keV, a = 0.1 keV), their equivalent widths are 



Wxa = 98^5 eV, Wk^ < 8eV, and Wcs = 22+^ eV, respectively. We note that, in this fit, the 



intrin sic widths of the iron lines were taken from their Chandra/HETG values (lYaqoob et al. 



20051 ). Substituting the individual Gaussian lines and pexrav component for the pexmon 
model results in a modest change in the global goodness-of-fit (Ax^/Az/ = —32/ — 3), likely 
owing to the free normalizations of the emission lines (they are set at fixed ratios within 
pexmon). No statistically significant change is seen in the other model parameters. 

For illustrative purposes only. Fig. |2t shows the residuals that remain when the broad 
iron line is removed from this spectral model, and the remaining model parameters are 
re-fit. An obvious broad line remains. However, there are two reasons why this cannot 
be interpreted as the "the broad line profile" for this object. Firstly, NGC 3783 has a 
well-known, high-column density warm absorber that, while principally affecting the soft 
spectrum, can also introduce subtle spectral curvature up to lOkeV or more. Secondly, the 
broadened iron line is just the tip of the iceberg; especially when the accretion disk is ionized, 
the rest of the reflection spectrum has a sub-dominant but significant contribution that must 
be considered. The statistically unlikely value of the disk inclination derived from the simple 
fit {i < 9°) is a signal of these issues. For these reasons, we are forced into global modeling 
of the full 0.7 — 45keV spectrum. 



•^The pexmon model ( Nandra et al.l 120071) is a modification of the commonly used pexrav model 
(jMagdziarz fc Zdziarskil 119951 ) which, in addition to the Compton backscattered reflection continuum, also 
models the Ka and K/3 emission lines of iron, the Compton shoulder of the iron Ka line, and the Ka line of 
nickel. The lines are included at the appropriate normalization for the assumed inclination, abundance, and 
reflection fraction. Hence, pexmon is superior to the usual "pexrav+gaussian" model since the strengths 
of the Compton reflection continuum and fluorescent lines are forced to be self-consistent. Assumptions do 
need to be made, however, when employing this model. In particular, we fix the inclination parameter of 
this component to be i = 60° and assume cosmic abundances. 
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3.2. Guidance from the long Chandra/HETG observation 



It is well known that NGC 3783 possesses a high-column density warm absorber (WA; 



e.g., iReynoldsl 119971 ): this is the greatest complexity we face when modeling the X-ray spec- 
trum of this source. For guidance, we turn to the long (900 ks) observation of NGC 3783 
with the HETGS on Chandra. Extensive analyses of the H ETG data have been published 
( iKaspi et al.ll2002l : iKrongold et al.ll2003l : iNetzer et al.ll2003l ): however, to retain consistency 
and utilize the latest calibrations, we have retrieved these data from the tgcat database! 
and have reanalyzed the Ist-order MEG+HEG spectra. 

In more detail, we obtain the ChandraRETG-daXa for NGC 3783 from tgcat for each 
of the OBSIDs corresponding to the 900 ks campaign and, coadd together spectra for a given 
order of a given grating. As a result, we obtain four spectral files corresponding to the time- 
average ±lst order spectrum from each of the HEG and the MEG. These were binned to a 
minimum of 15 photons per spectral bin in order to validate the use of techniques while 
still maintaining spectral resolution. We then jointly analyzed these spectra, noticing the 
0.5-7 keV range in the MEG data and the 1-7.5 keV range in the HEG data. We permitted 
the overall cross-normalization between these four spectra to be free parameters; in all cases, 
the best-fitting cross-normalization is within 5% of unity. 

Fitting these data with a power-law modified by the effects of Galactic absorption 
[Nh = 9.91 X 10^°cm~^ ; described using the phabs model of XSPEC) results in a very poor 
fit with x^/^ = 58962/13112(4.50). The residuals suggest a soft excess component, soft 
X-ray absorption by a WA, and a prominent fl uorescent iron Kg line a 6 .4 keV. The effect 
of the WA is modeled using the XSTAR code (IKallman &: Bautistall200l[ ): for an absorber 
of a given column density A'^^^ and ionization parameter ^, XSTAR is used to compute the 
absorption imprinted on a power-law X-ray spectrum. We compute a grid of XSTAR models, 
logarithmically sampling a range of column densities in the range Nh '■ 10^^ — 10^^ cm~^ and 
a range of ionization parameters in the range ^ : 1 — 10^ erg cm s~^, for use in spectral 
fitting. In the construction of the WA grids it is assumed that elemental abundances are 
fixed to solar valued, and that the turbulent velocity of the WA is 200kms~^ . Dramatic 
improvements in the goodness-of-fit are found by the inclusion in the model of three zones 
of WA. To begin with, each WA component is included assuming that the absorbing gas is 
at rest with respect to NGC 3783; the improvement in the fit upon the addition of each WA 
component was Ax^ = —26316, —3096 and —2490. The inclusion of a fourth zone led to a 



''http://tgcat.niit.edu/ 

^ http: / /heasarc. nasa.gov/lheasoft / xstar/docs/html / node41 .html 
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much smaller improvement in the fit and hence was deemed inappropriate. The residuals 



from the 3-WA fit do indicate a soft excess. Following iKrongold et al.l ( 120031 ). we model the 
soft excess with a blackbody component (this is intended to be a phenomenological, not a 
physical, description of the soft excess; see discussion in §3.3p resulting in an improvement 
in the fit of A^VAz/ = -480/ - 2 (i.e., x^l^ = 26582/13104 (2.03)). 

While providing a decent fit to the global spectrum, the model thus described leaves 
prominent unmodeled emission and absorption lines, the most prominent of which is the 
iron fluorescent emission line at 6.4 keV. Fitting the iron line with a simple Gaussian model 
improves the goodness-of-fit by Ax^/Ai/ = —593/ — 3, with a line energy E = 6.398 ± 
0.002 keV (confirming the identification of cold iron-Ka), FWHM= 2000 ± 300 km s"^ and 
equivalent width Wko = 88±6eV. However, since we believe that this component originates 
from reflection, we shall henceforth model it using pexmon; replacing the simple Gaussian 
with the pexmon model convolved with a Gaussian velocity profile (with FWHM= 1800 ± 
300 km s^^ ) results in a slightly better fit (A^^ = —18). At the soft end of the spectrum, 
the Ka emission triplet of Ovil (at 0.574 keV, 0.569 keV, 0.561 keV) as well as the Ka 
emission line of O vill (at 0.654 keV) are clearly visible. Modeling these as Gaussian lines at 
the redshift of NGC 3783 with common velocity width yields a further improvement in the 
goodness of fit (Ax^ = —218) with best-fitting FWHM= 700 ± 150kms~^ and equivalent 
widths 1^0.574 = 26 ± 6 eV, 1^0.569 = 14 ± 6 eV, 1^0.561 = 47 ± 9 eV, and 1^0.654 = 23 ± 5 eV. 

It is well known that the WA in this and many other objects corresponds to outflowing 
gas. Relaxing the constraint that the WA zones are at the systemic redshift of NGC 3783 
yields a large improvement in the fit (AxV'^'^ = 4247/ - 3; xV'^ = 21532/13094(1.64)), 
with implied line-of-sight outflow velocities in the 500 — 1000 km s^^ range. These velocities 
as well as the other parameters defining the best-fit model for the HETG data are listed 
in Table 1. The spectral model described in this section (power-law continuum, three- zone 
WA, blackbody soft excess, reflection from distant neutral material, and emission lines from 
O VII and O vill) describes the vast majority of spectral features seen in the HETG data 
(see Fig. [3]). 



3.3. Global Modeling of the 0.7-45 keV Spectrum 

To extract the maximal information from the full-band (0.7-45 keV) Suzaku spectrum 
of NGC 3783, we must compare the data to a global spectral model which is as physically 
self-consistent and realistic as possible. In constructing this global model, we draw guidance 
from our heuristic analysis of the hard-band spectrum ( §3.ip as well as the results from 
the C/ian(/ra/HETG ( §3.2p . The primary continuum emission is taken to be a power-law 
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Fig. 3. — Folded Chandra/IIETG spectrum and best fitting model as a function of observed energy. 
As described in the text ( §3.2p . the model is fitted simultaneously to the ±1 MEG (0.5-7 keV) and 
HEG (1-7.5 keV) data. However, for clarity, we only show here the — 1-order MEG data (black; 
first three panels) and the —1-order HEG data (blue; bottom panel). 
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5 10 20 5 10 

Energy (keV) Energy (keV) 

Fig. 4. — Results of fitting tlie XIS+PIN spectrum with a model that includes the warm absorbers, 



distant reflection and scattering/leaked soft component but not the relativistic ionized accretion 
disk. While the fitting is performed on the 0.7-45 keV spectrum, we show for clarity only the 
residuals above 3 keV. Left: Strong residuals indicative of a broad iron line and Compton reflection 
hump are clearly visible. This motivates the inclusion of a relativistic disk component into the 
spectral model. The XIS 0+3 data are shown in black and XIS 1 data are in red, while the 
HXD/PIN data are in blue. The solid green line represents a data-to- model ratio of unity. Right: 
Zoom-in on the Fe K line region. 



(photon index F) with a soft excess which we describe as a blackbody (temperature T). 
X-ray reflection of this continuum from cold, distant material (possibly associated with the 
dusty/molecular torus of unified Seyfert schemes) is described using the pexmon model (see 
§3.ip . As discussed in §3.H the inclinatio n of the pexmon is fi xed at z = 60° and the abun- 
dances are fixed to be solar, as defined in iNandra et al.l ( 2007 ). These emission components 
are then absorbed by a three-zone WA modeled using the XSTAR tables described in §3.2( 
the column density N\y and ionization parameter ^ of each zone is taken to be a free pa- 
rameter rather than being fixed to the HETG value. Since the Suzaku/XlS detectors do not 
have the spectral resolution capable of constraining the outfiow velocities of the various WA 
zones, we have elected to hold the redshifts of these components fixed at the cosmological 
value for NGC 3783. Statistically indistinguishable results are obtained if we, instead, fix 
the outfiow velocities to the HETG-derived values. For completeness, we also allow for some 
fraction /sc of the continuum to be scattered around (or leak through) the WAs, i.e., our 
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Energy (keV) Energy (keV) 

Fig. 5. — Global modeling of the 0.7 - 45keV XIS-FI+PIN data. The left panel shows the 
resulting residuals from fitting the model (including the relativistic accretion disk) discussed in 
^3.31 Data point colors are as in Fig. 2. The right panel shows the best-fitting model color coded 
as follows: (a) green line, continuum power-law emission; (b) dark blue line, cold and ionized iron 
line emission from distant matter; (c) red line, soft excess modeled as blackbody; (d) magenta line, 
significant emission that scatters around or leaks through the warm absorber; (e) light blue line, 
relativistically-smeared disk reflection, (f) thick black line, total summed model spectrum. Warm 
absorption affects all components apart from (d). 



model allows for "partial covering" . 



Fitting this model to the 0.7-45 keV Suzaku data results in a poor fit (x^/ v = 1206/679 (1.78)) 
and strong residuals which indicate the presence of the broad iron line as well as addi- 
tional reflection beyond that associated with the narrow iron line (Fig. |4]). This leads us 
to include relativistically smeared reflection from an ionized accretion disk into the spectra. 1 



model; operationally, we use the ionized reflection model reflionx (IRoss fc Fabian 



convolved with the variable spin relativistic smearing model relc onv (Dauser et al 



20051) 

2om 



The r elconv model is a further evolution of the kerrconv model of iBrenneman fc Reynolds 
( 120061 ) ■ employing faster and more accurate line- integration schemes and allowing black 
hole spin to be fit as a free parameter for prograde, non-spinning and retrograde spins 
(a e [—0.998,0.998]). While the fit achieved with this blurred reflection model is not sta- 
tistically ideal (x^/i^ = 917/664(1.38)), there are no broad-band residuals (Fig. and 
much of the contribution to the excess originates from fine details of the WA-dominated 
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region below 1.5 keV. The model is shown in Fig. |5]3, and the best-fitting parameter values 
are shown in Table 1. 

The parameters defining the best-fitting model for the 0.7 — 45keV data are shown in 
Table [1] along with their 90% confidence ranges. Under the assumption that we can identify 
the low-, medium-, and high-ionization components seen in the 2001-HETG observation 
with those seen in out 2009- Suzaku data, we see that both the column density and the 
ionization state of the low-ionization absorber have increased somewhat (Alog^ ~ 0.28, 
A^WA ~ 3 X lO^^cm"^). In contrast, the medium- and high-ionization absorbers have 
slightly dropped in ionization parameter. Given that these different zones are likely at very 
different distances from the central engine with very different plasma densities, they will 
possess very different recombination/photoionization timescales and hence will respond to 
changes in the ionization flux on different timescales. Thus, it is not surprising that we see 
a mixture of increasing and decreasing ionization states in the various WA zones. 

We also note a change in the both the temperature and normalization of the blackbody 
component between the 2001-HETG and 2009-Suzaku data. This merits some discussion. 
The blackbody co mponent used t o phen omenologically parameterize the soft excess was 
flrst employed by iKrongold et al.l (120031 ) . who found kT = 0.10 ± 0.03 keV and Abb = 



2.0 ± 0.7 X 10 ^, where the normalization is in units of L^i^/ D\q (L39 is luminosity of the 
component in units of lO^^ergs"^ and Diq is distance to the source in units of lOkpc). Our 



analysis of the same HETG data conflrms the lKrongold et al.l (120031 ) result. By contrast, our 
analysis of the S'uza/En/XIS+PIN spectra finds a lower temperature {kT = 0.060!l^4 keV) and 
a normalization that is almost two orders of mag nitude greater (Abb = 8.4^^;? x lO'^). We 
stress that neither the use of a blackbody to model the soft excess nor the precise change in 
the parameters of the blackbody should be interpreted literally. In particular, the significant 
change in the normalization of this component is misleading — the lower energy cutoffs in 
both the HETG analysis (0.5 keV) and the XIS analysis (0.7keV) are much higher than the 
peak of this blackbody component and, thus, only the Wien tail of this component is playing 
any role in the spectral fitting. Given this fact, even a modest drop in the temperature must 
be compensated for by a large increase in normalization in order to have a comparable 
contribution in the observed energy band. While the physical nature of the soft excess is 
of intrinsic interest, it is beyond the scope of this paper. We have verified that different 
treatments of the soft excess (replacing the blackbody spectrum with bremsstrahlung or a 
steep power-law component) do not affect the interpretation of the spectrum above 2 keV. 

The principal focus of this work is the signature of the relativistic accretion disk. Our 
global fit finds refiection from a rather low-ionization accretion disk (^ < 9 erg cm s~^ ) ex- 
tending down to the innermost stable circular orbit (ISCO) of a rapidly rotating black hole 



- 14 - 



(a > 0.98). The emissivity/irradiation profile defining the reflection spectrum, modeled as 
a broken power-law, is found to have an inner power- law index of qi = 5.2lQg breaking to 
q2 = 2.9 ± 0.2 at a radius of = 5Atl'^rg. If these indices are tied together in the model, 
i.e., if qi = q2, the fit worsens considerably (Ax^/Az/ = +21/ + 1, with qi = q2 = 3.0 ± 0.3) 
and the black hole spin is also less tightly constrained: a > 0.25. 

The iron abundance of the disk has been constrained to lie between 2.8 — 4.6 times 
solar. To probe the robustness of this constraint we have refitted the data in three different 
ways, each allowing for slight differences in the way the iron abundance was handled in 
our model: (1) fixing Fe/solar of the distant reflector (pexmon) to that of the inner disk 
reflection (reflionx), with both values frozen at Fe/solar= 1; (2) allowing these linked 
abundances allowed to vary freely; and (3) allowing both abundances to vary freely and 
independently. Compared with the global best-fit, scenario (1) resulted in a worsening of 
goodness-of-fit by Ax^ = 33 and unconstrained black hole spin at the 90% confidence level, 
scenario (2) resulted in a marginal decrease in goodness-of-fit by Ax^ = 7 (no change in 
spin constraints) and scenario (3) yielded no change in the goodness-of-fit (no change in spin 
constraints). In summary, the high iron abundance of the reflionx is statistically preferred 
to the solar value; the high spin value is dependent upon the high iron abundance, but the 
high abundance is strongly preferred in the fit. Because the iron abundance of the distant 
reflector could not be constrained independently of the relativistic reflector, the pexmon and 
reflionx iron abundances have been linked in our best model fit. 

To gauge the importance of the (complex) soft spectrum on our global fit, we have also 
conducted a restricted hard-band (3-45 keV) fit. Since a hard-band fit cannot constrain the 
parameters of the WA or soft excess, these parameters are constrained to lie within their 90% 
confidence ranges as derived from the 0.7-45 keV analysis. To be most conservative, we also 
relax the constraint on the XIS/PIN cross-normalization, allowing it to be a free parameter. 
The resulting fit is listed in the last column of Tabled! For this fit x^/z/ = 499/527(0.95), 
a great improvement over the 0.7 — 45keV fit, and confirmation that the small residuals 
below ~ 1.5 keV are the primary contribution to the large reduced oi the full spectral fit. 
While the parameter values are equivalent to those of the 0.7 — 45keV fit within errors, the 
uncertainties on the parameters are larger when only the hard spectrum is considered. This 
is especially true for the inner disk emissivity and break radius of the relconv model, which 
exhibit a strong degeneracy without the soft spectrum data. Figure [6] shows the confidence 
contours on the (gi, g2)-plane for this hard band fit; we see that the outer disk emissivity 
index, q2, is well-constrained, whereas the constraints on the inner disk emissivity index, 
gi, are clearly worse. Fixing the XIS/PIN cross- normalization at the nominal value of 1.18 
tightens the constraints but still leaves a significant degeneracy between qi and rbr. (Fig. [6]). 
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Using the best-fitting 0.7 — 45keV spectral model with the XIS normalization, the 
2 - lOkeV observed-frame flux of NGC 3783 is F2-10 = 6.04 x 10^^^ ergcm'^ s"^ . Adopting 
a standard cosmological model {Hq = 71kms~^ Mpc^^, Hm = 0.3, = 0.7), this implies a 
rest-frame luminosity of -^2-10 = 1-26 x 10^'^ erg s~^ . The hard X-ray band yields a 16—45 keV 
flux of -F16-45 = 1-07 X 10~^° erg cm~^ s^^ for a rest-frame luminosity of -L16-45 = 2.24 x 
10"^^ ergs"^ . 

3.4. The Spin of the Black Hole 

Our flducial spectral model discussed above yields a spin constraint of a > 0.98 (90% 
confldence), or a > 0.88 (99% confldence). However, given the subtle nature of the spin 
measurements, it is useful to address the systematic issues that may be introduced by the 
modehng and analysis techniques. 

We can assess the role of different analysis-related assumptions on the derived spin by 
comparing the variation of with a. Figure [7] (black line) shows Ax^(a) = x^(cj) ~ Xbest-fit 
from our flducial analysis that underlies the constraint just quoted. It is interesting to note 
the non-monotonic nature of the x^-space above a ~ 0.75. We consider a few variants from 
this flducial analysis in order to probe the sensitivity of the spin measurement. An important 
issue is the extent to which the warm absorber parameters are trading-off with the derived 
black hole spin. Thus, we repeat the analysis with the warm absorber parameters (column 
densities and ionization parameters for all three zones) flxed at their best-flt values from the 
flducial model. The resulting spin constraints are shown in Fig. [7] (a > 0.98; red line) and 
are very similar to the flducial model, indicating little or no degeneracy between spin and 
the warm absorber parameters. 

Secondly, to the extent that the strength of the Compton reflection hump is important, 
we may be concerned about the effect of cross-calibration errors in the flux normalization 
between the XIS and the PIN spectra. Thus, we repeat the spectral analysis, leaving the 
cross-normalization factor as a free parameter. The best-flt value is slightly smaller than our 
flducial value (1.15 vs. 1.18), but the improvement in the goodness-of-flt is only marginally 
signiflcant (Ax^ = 6 for one additional degree of freedom). The spin is constrained to be 
slightly smaller than that of the flducial model (to 90% confldence, a = 0.92 — 0.95; Fig. [TJ 
blue hne). 

Lastly, we may be concerned that the spin flts are being driven by the contribution of 
the ionized disk to the high-S/N by highly-complex region of the spectrum below 1.5 keV. 
Thus, we have repeated our analysis including data only above 3keV. Here, too, we allow 
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Table 1: Spectral fit parameters. All errors are quoted at the 90% confidence level for one interesting 
parameter (Ax^ = 2.7). Parameters marked with an "(f)" had their values fixed during the fit. 
Units of normalization are in phcm-^s"^, column density is in units of lO^'^ cm-^ , ionization 
parameter is in erg cm s-^ , iron abundance is relative to solar (linked between the reflionx and 
pexmon reflection components), inclination is in degrees, radii are in rg, and spin parameter is 
dimensionless, but is defined as a = cJ/GM^. See §3 for details. 
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the XIS/PIN cross-normalization to be a free parameter. Given the lack of data at soft 
energies to constrain them, the WA, BB and scattered fraction components were frozen to 
their best-fitting values for the full-band, free cross-normalization case (which is, within 
errors, identical to the WA parameters for the fiducial model). Yet again, the best-fit spin 
parameter is similar to that of the fiducial model (a > 0.95; Fig. [TJ green line). This indicates 
that the fitted spin value is indeed driven by the Fe K band. 



4. Discussion and Conclusions 

The X-ray spectrum of NGC 3783 is complicated; in addition to the effects of a multi- 
zone warm absorber, there are suggestions that some fraction (17%) of the primary X- 
ray emission can scatter around or leak through the warm absorber. However, despite 
this complexity, the high-S/N and broad bandpass of Suzaku allows us to robustly detect 
and study the relativistically-smeared X-ray reflection spectrum from the surface of the 
inner accretion disk. Assuming that the region within the ge neral relativistic radius o f 



marginal stability does not contribute to the reflection spectrum ([Reynolds fc Fabian! l2008l ) 
we determine a lower limit of a > 0.98 (90% confidence) to the dimensionless spin parameter 
of the black hole. Even at the 99% confidence level, we can constrain the spin to be a > 0.88. 
Relaxing the assumed XIS/PIN cross-normalization or neglecting the soft-band data (but 
then freezing the WA parameters) allows the model to find a slightly better fit and makes 
the constraints slightly lower (a = 0.92 — 0.95, a > 0.95 at 90% confidence, respectively; 
a > 0.88, a > 0.90 at 99% confidence, respectively). 

Including this result, four out of eight of the AGN with reliable spin measurements 
may have spins greater than a = 0.8 (see Table 2). Spin measurements for more sources 
are required before we can draw any conclusions about the spin distribution function, but 
here we note that there are potentially important selection effects biasing any flux-limited 
sample towards high spin values. For standard accretion models, the efficiency of black 
hole accretion increases as the spin of the black hole increases. So, all else being equal, 
an accreting, rapidly spinning black hole will be more luminous than an accreting, slowing 
spinning black hole and hence will be over-represented in flux-limited samples. 

We illustrate this effect by calculating the selection bias given some very simple as- 
sumptions. Suppose that a flux limited sample is constructed in some band B. The accretion 
luminosity in that band will be given by 

L = KbvMc^ (1) 



where Kb is the fraction of luminosity appearing in band B (i.e. the reciprocal of the 
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Table 2: Summary of black hole spin measurements derived from relativistic reflection fitting of SMBH spectra. Data are taken 
with Suzaku except for 1H0707-495, which was observed with XMM-Newton, and MCG-6-30-15, in which the data from XMM 
and Suzaku are consistent with each other. Spin (a) is dimensionless, as defined previously. Wxa denotes the equivalent width 
of the broad iron line relative to the continuum in units of eV. Parameter qi represents the inner disk emissivity index and is 
unitless. Fe/solar is the iron abundance of the inner disk in solar units, while ^ is its ionization parameter in units of ergcms"^ . 
M is the mass of the black hole in solar masses, and -Lboi/-^^Edd is the Eddington ratio of its luminous output. Host denotes the 
galaxy host type and WA denotes the presence /absence of a warm absorber. Values marked with an asterisk either were fixed 
in the fit or have unknown errors. Al l masses are frq r n iPeterson et al.l (120041) exc ept MCG-6-30-15, 1H070 7-495 and SWIFT 
J2127.4+5654, which are taken from iMcHardy et al.l (|2005l ). Izoghbi et all (|201ol ) and iMahzia et al.l (|2008l ). respectively Ah 
bolometric luminosities are from lWoo &: Urrvl (j2002l ) except for the same three sources. The same references for MCG-6-30-15 
and SWIFT J2127.4+5654 are used, but host types for 1H0707-495 and SWIFT J2127.4+5654 are unknown. 



'^Brenneman & Reynolds (2006), Miniutti et al. (2007). 

''Schmoll et al. (2009), though note some discrepancies with Patrick et al. (2010). 
'^Miniutti et al. (2009), though note some discrepancies with Patrick et al. (2010). 
■^Zoghbi et al. (2010), de La Calle Perez et al. (2010). 
^Gallo et al. (2005, 2010). 
^Patrick et al. (2010). 
^This work. 
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bolometric correction), rj is the accretion efficiency, and M is the mass accretion rate. Now 
let us assume that M has no exphcit spin dependence (e.g. is determined by the larger 
circumnuclear environment), and that the spectral energy distribution and hence Kb is 
independent of spin. Thus, the space density of sources with accretion rates in the range 
M — )■ M + dM and spins in the range a — )■ a + da, denoted $(M, a) dM da, can be taken as 
a given function set by the astrophysics of black hole growth. 

We assume a Euclidean universe, valid for the local/bright AGN samples relevant for spin 
measurements with the current generation of X-ray observatories. The number of sources 
in a flux-limited sample with luminosity in the range L ^ L + dL and spins in the range 
a ^ a + da is then, 

dN o^^{L,a)L^/^ dLda, (2) 

where a) dL da is the space density of sources with luminosity in the range L ^ L + dL 
and spins in the range a — > a + da. Transforming into the (M, a)-plane gives, 

dN (x^{M,a)7f''^dMda. (3) 

Using our assumption that the mass accretion rate is independent of spin, we can separate 
$(M , a) into an accretion rate dependent space density n{M) and a spin distribution function 
/(a), $(M, a) = n{M)f{a). We can then integrate Eqn. [3]over M in order to determine the 
number of sources in a fiux-limited sample with spins in the range a — )■ a + da: 

dNa oc /(a)?7(a)2/2 da (^j n{M) dM^ . (4) 

For illustration purposes, let us examine eqn. H] in the case of a completely fiat spin distribu- 
tion where /(a) = constant for a G [0,amax] and is zero otherwise. Thus, half of the parent 
population as a whole has a > a^^^/2. We find that if amax = 0.95, then half of sources in 
the fiux limited sample will have a > 0.67; for amax = 0.99 we find that half of the sources 
in the sample have a > 0.73. 

Generalizing away from a fiat spin distribution, we can consider spin distribution func- 
tions of the form /(a) oc a^. Within this simple framework, we require /(a) oc a (i.e. p = 1.0) 
in order to produce fiux-limited samples where half of sources have a > 0.84 (assuming 
Omax = 0.95). For high-spin weighted distribution functions such as this, the selection bias 
is stronger; only 20% of objects in the volume-limited parent sample actually have a > 0.84. 
Of course, given the small number statistics and highly inhomogeneous selection functions 
for the current spin measurements, it is too early to draw any conclusions about the need 
for a high-spin biased distribution function. 

We are extremely grateful to our NASA and J AX A colleagues in the Suzaku project for 
enabling these Key Project data to be collected. We thank Martin Elvis and Cole Miller for 
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Fig. 6. — Confidence contours on the {qi,q2)-plsine for tlie 3-45 keV band fit assuming a XIS/PIN 
cross-normalization that is free (left panel) or fixed at the nominal value of 1.18 (right panel). Blue, 
red and green lines show the 68%, 90% and 99% confidence contours for two interesting parameters 
(Ax^ = 2.3, 4.6, 9.2), respectively. The dot-dashed black line represents qi = q2- 
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Fig. 7. — Goodness-of-fit parameter Ax^ as a function of the assumed black hole spin for our 
fiducial model (black line), the fiducial model except with frozen WA parameters (red line), the 
fiducial model except with free XIS/PIN normalization (blue line), and a hard-band (> 3keV) fit 
only (green line). Confidence levels are indicated with horizontal black lines and are derived for 
one interesting parameter. See text in §3.4 for details. 



